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Abstract 

This  paper  presents  the  results  of  the  modeling  and  analysis  of  the  catalyst  layer  of  a  PEM  fuel  cell  cathode.  We  nondimensionalized 
the  governing  equations  and  introduced  three  dimensionless  parameters,  tt\,  jtk  and  ny,,  to  simplify  the  description  of  the  complicated 
phenomena  occurring  in  the  catalyst  layer.  The  three  dimensionless  parameters  indicate  the  resistance  to  oxygen  reduction,  the  resistance 
to  proton  conduction,  and  the  resistance  to  oxygen  diffusion,  respectively.  The  profiles  of  oxygen  concentration,  reaction  rate,  and  current 
density  in  the  catalyst  layer  were  investigated  in  terms  of  jtk  and  ny,.  It  was  found  that  the  dimensionless  overpotential  across  the  catalyst 
layertakestheformof  i)cat  =  —In  n\— /(jtk,  ny, ).  The  functional  relation  of  /(jtk,  ny,),  which  represents  the  mixture  of  ohmic  overpotential 
and  concentration  overpotential,  was  then  evaluated  for  various  combinations  of  jtk  and  jtd,  and  finally  two  correlation  equations  for  i) cat 
were  proposed.  The  plot  of  f{n k,  ny,)  and  the  proposed  correlation  equations  are  helpful  for  analysis,  prediction  and  optimization  of  the 
performance  of  a  PEM  fuel  cell. 
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1.  Introduction 

The  cathodic  overpotential,  which  constitutes  the  largest 
irreversible  losses  in  the  cell  voltage,  is  influenced  by  sev¬ 
eral  physical  and  operating  parameters  such  as  the  cell 
current  density,  the  active  catalyst  surface  area,  conduc¬ 
tivities  and  thickness  of  the  catalyst  layer,  the  concen¬ 
tration  and  diffusion  coefficient  of  oxygen,  and  so  on. 
Mathematical  models  are  a  useful  tool  for  analysis,  pre¬ 
diction  and  optimization  of  the  performance  of  the  fuel 
cell  cathode.  Previous  studies  have  offered  many  mod¬ 
els  for  investigating  the  phenomena  and  overpotential 
occurring  in  the  catalyst  layer  of  a  fuel  cell  electrode. 
The  mathematical  models  devoted  to  the  catalyst  layer 
of  a  fuel  cell  electrode  include  the  pseudo-homogeneous 
model  [1-10],  the  single  pore  model  [11],  the  thin  film 
model  [12],  and  the  agglomerate  model  [4,12].  They  in¬ 
corporate  different  degrees  of  complexity  and  applicability. 
Among  these  models,  the  pseudo-homogeneous  model  is 
the  one  most  often  used,  mainly  because  it  is  straightfor¬ 
ward  and  its  outcome  compares  well  with  experimental 
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results.  The  pseudo-homogeneous  model  assumes  that 
the  catalyst  layer  has  a  uniform  thickness  and  uniform 
distribution  of  its  various  components,  i.e.  the  sizes  of 
catalyst  particles,  ionomer,  and  gas  pores  are  extremely 
small  and  are  uniformly  distributed  throughout  the  catalyst 
layer. 

The  pseudo-homogeneous  model  has  been  adopted  by 
the  following  studies  in  regard  to  the  PEM  fuel  cell  cathode. 
Bernardi  and  Verbruge  [1,2]  presented  a  complete  model 
for  a  gas  diffusion  cathode  bonded  to  a  polymer  electrolyte. 
They  incorporated  the  electro-osmotic  and  pressure-driven 
water  transport  consonant  of  the  PEM  fuel  cell  to  investi¬ 
gate  the  overpotential  characteristics,  water  transport,  and 
catalyst  utilization.  Springer  et  al.  [3]  developed  a  model 
considering,  in  detail,  the  voltage  losses  caused  by  inter¬ 
facial  kinetics  at  the  Pt/ionomer  interface,  gas  transport 
and  proton  conduction  limitations,  and  gas  transport  lim¬ 
itations  in  the  cathode  backing.  Broka  and  Ekdunge  [4] 
studied  the  discrepancies  in  polarization  curves  given  by 
the  pseudo-homogeneous  model  and  the  agglomerate  model 
under  various  values  of  oxygen  permeability,  conductivities 
and  thickness  of  the  catalyst  layer.  Marr  and  Li  [7]  used  a 
mathematical  model  to  study  the  performance  of  the  cath¬ 
ode  catalyst  layer  in  a  PEM  fuel  cell  by  including  both 
electrochemical  reaction,  and  mass  transport  processes. 
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The  composition  and  performance  optimization  of  catalyst 
platinum  has  also  been  investigated. 

All  previous  studies  on  the  pseudo-homogeneous  model 
for  the  catalyst  layer  of  a  PEM  fuel  cell  cathode  have  taken 
into  account  all  the  parameters  influencing  the  performance 
of  a  fuel  cell,  and  investigated  their  influences  separately. 
Owing  to  the  fact  that  the  number  of  influential  parame¬ 
ters  is  large,  the  analysis  and  optimization  of  a  PEM  fuel 
cell  cathode  turn  out  to  be  very  complicated.  However,  the 
influential  parameters  can  be  organized  into  a  lesser  num¬ 
ber  of  dimensionless  groups,  making  the  problem  easier  to 
deal  with.  In  the  present  study  we  nondimensionalized  the 
governing  equations  and  boundary  conditions,  to  reduce  the 
description  of  the  complicated  phenomena  occurring  in  the 
catalyst  layer  of  a  PEM  fuel  cell  cathode  to  three  dimension¬ 
less  parameters.  The  cathodic  overpotential  was  evaluated 
for  various  combinations  of  the  dimensionless  parameters, 
and  then  the  relationship  between  the  cathodic  overpotential 
and  the  dimensionless  parameters  was  obtained. 


2.  Model  description 

Fig.  1  presents  a  schematic  of  the  cathode  bonded  to  a 
proton  exchange  membrane  in  a  PEM  fuel  cell.  The  system 
can  be  divided  into  three  regions:  a  membrane  region  of 
solid  polymer  electrolyte  (PEM),  an  active  catalyst  region 
(catalyst  layer)  that  provides  a  catalytic  site  for  the  reduc¬ 
tion  of  oxygen,  and  a  diffusion  region  (gas  diffusion  layer) 
composed  of  highly  porous  and  conductive  material. 

The  gas  diffusion  layer  is  located  adjacent  to  the  air/C>2 
flow  channel  and  the  current  collector.  Its  open  pores  and 
electronically  conductive  material  transport  O2  and  electrons 
to  the  catalyst  layer,  respectively.  The  PEM  is  an  electronic 
insulator,  but  an  excellent  conductor  of  protons.  The  pro¬ 
tons  produced  at  the  anode  go  through  the  PEM  to  the  cat¬ 
alyst  layer  of  the  cathode,  where  they  react  with  electrons 
and  O2  to  form  H2O.  The  catalyst  layer  is  formed  as  a  thin 
layer  of  proton  conductive  ionomer  (e.g.  Nafion®),  and  a 
carbon-supported  catalyst  (e.g.  Pt/C).  The  ionomer  portion 
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Fig.  1.  Schematic  of  the  cathode  bonded  to  a  proton  exchange  membrane 
in  a  PEM  fuel  cell. 


of  the  catalyst  layer  forms  the  transport  avenue  for  protons, 
and  the  matrix  portion  consisting  of  carbon  and  catalyst 
plays  the  role  of  electronic  conductor.  Studies  on  the  mi¬ 
crostructure  of  the  catalyst  layer  have  shown  the  existence 
of  void  space  forming  a  network  of  gas  channels  that  dom¬ 
inate  the  transport  of  O2  to  the  active  site  [3,4,6]. 

In  the  present  study  we  focus  on  the  catalyst  layer  of  the 
cathode,  which  is  in  the  region  of  0  <  z  <  <5C.  In  model¬ 
ing  this  catalyst  layer,  the  following  assumptions  have  been 
adopted. 

•  The  fuel  cell  is  operating  under  a  steady  state. 

•  The  catalyst  layer  is  isothermal. 

•  The  distribution  of  oxygen  concentration,  current  density, 
and  the  overpotential  in  the  catalyst  layer  is  considered  to 
be  one  dimensional. 

•  The  oxygen  permeation  through  the  PEM  is  negligible. 

•  The  catalyst  layer  can  be  regarded  as  a  pseudo-homoge¬ 
neous  film,  so  that  the  active  surface  area  of  catalysts, 
void  space  and  ionomer  are  uniformly  distributed  in  it. 

•  The  void  space  is  sufficiently  large  so  that  the  Knudsen 
diffusion  is  negligible  [8,13].  Besides,  the  convective  mass 
transfer  of  oxygen,  driven  by  liquid  water  and  water  vapor, 
is  unimportant  compared  to  the  bulk  diffusion. 

•  The  electronic  conductivity  of  Pt/C  is  substantially  greater 
than  the  protonic  conductivity  of  ionomer  [1,3, 6, 7], 


The  rate  of  the  electrochemical  reaction  in  the  catalyst  layer 
can  be  described  using  the  Butler- Volmer  rate  expression. 
This  is  then  simplified  to  give  a  Tafel-type  equation  in  terms 
of  the  oxygen  concentration  as 


d  i  .  Co, 

—  =  AviQ7^-  exp 


dz 


r 


02 


/  acFr]\ 

V  RT  ) 


(1) 


In  Eq.  (1),  i  is  the  local  protonic  current  density,  Av  the 
specific  area  of  the  active  surface,  i'o  the  reference  exchange 
current  density,  Co,  the  local  oxygen  concentration,  C^h  the 
reference  oxygen  concentration,  which  is  associated  with 
/(),  ac  the  cathodic  transfer  coefficient,  and  11  is  the  local 
cathodic  overpotential. The  local  cathodic  overpotential  r](z) 
is  related  to  the  potential  difference  of  0s(z)  and  0m(z)  by 
the  equation: 

n(z)  =  [0s  (z)  -  0m  (z)]  -  0ref  (2) 


where  0S  is  the  potential  of  the  electronically  conductive 
phase  in  the  catalyst  layer  (i.e.  the  solid  portion  composed 
of  carbon  and  catalyst  particles),  0m  the  potential  of  the 
ionomer  phase,  and  0ref  the  reference  potential  difference. 
Both  0s  and  0m  increase  in  the  positive  z-direction,  and 
equations  of  Ohm’s  law  for  each  phase  can  be  expressed  as 
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where  I  is  the  operating  cell  current  density,  K*f f  and 
denote  the  effective  conductivities  of  the  solid  portion  and 
the  ionomer  phase,  respectively.  Therefore,  the  gradient  of 
the  overpotential  in  the  catalyst  layer  is  given  by 

dij  _  i  /I  1  \  . 

dz  ~~  Kff1-  +  Kf)1  (5) 


Because  is  much 
reduced  to 
d  ii  1 


greater  than  K^,  Eq.  (5)  can  be 

(6) 


dimensionless  parameters,  tti,  7Tk,  and  7Td,  the  governing 
equations  take  the  following  form: 


d  i 
dz 

d>] 
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(12) 


d  C 

=  -7TD(1  -  0 
dz 


(13) 


The  corresponding  boundary  conditions  become 


Note  that  the  cathodic  overpotential  rj(z)  and  its  first  deriva¬ 
tive  d;;/dz  are  both  negative.  Therefore,  r/(z)  is  a  negative 
and  monotonically  decreasing  function. 

One  can  relate  the  oxygen  concentration  with  the  protonic 
current  density  i  in  the  catalyst  layer  by  material  balance. 
On  the  assumption  that  both  convective  mass  transfer  and 
Knudsen  diffusion  of  oxygen  are  unimportant  when  com¬ 
pared  to  bulk  diffusion,  Fick’s  law  completely  describes  the 
mass  flux  of  oxygen  in  the  catalyst  layer.  Moreover,  as  oxy¬ 
gen  permeation  through  the  PEM  is  negligible,  the  mass  flux 
of  oxygen  crossing  any  location  z*  must  be  depleted  in  the 
region  z*  <  z  <  <5C.  Hence  the  material  balance  for  oxygen 
gives 


dCo2  I  —  i 

dz  ~~  4f7T[r 
02 


(7) 


where  /J^1  is  the  effective  diffusion  coefficient  of  oxygen 
in  the  catalyst  layer.  Both  and  l)f  can  be  calculated 
from  their  bulk  properties  and  volume  fractions  using  the 
Bmggeman’s  correction  [8,14,15],  i.e.  =  Km e^2 

and  =  Do2e3(~,  where  sm  is  the  volume  fraction 

of  ionomer,  and  s  is  the  volume  fraction  of  the  open 
pores. 

There  are  three  dependent  variables,  i,  Co,,  and  ?/,  in¬ 
volved  in  three  first-order  differential  equations,  Eqs.  (1),  (6) 
and  (7).  Adopting  the  oxygen  concentration  at  the  gas  dif¬ 
fusion  layer/catalyst  layer  interface  as  the  reference  oxy¬ 
gen  concentration,  the  corresponding  boundary  conditions 
are 


i(Sc )  =  I  (8) 

HO)  =  0  (9) 


Cq2  (0)  =  Cq2 


(10) 


3.  Solution  technique 

The  governing  equations  and  the  corresponding  boundary 
conditions  together  pose  a  boundary-value  problem.  This 
boundary-value  problem  can  be  reduced  using  the  following 
relations:  z  =  z/5c,  =  i/I,  C  —  Co2/Cq v  fj  =  rj/b, 

where  b  —  RT/oicF  is  the  Tafel-slope.  By  introducing  three 


k  i)  =  1  (i4) 

m  =  0  (i5) 


C(0)  =  1 


(16) 


The  three  it’s  are  meaningful  dimensionless  parameters  in 
the  characterization  of  the  catalyst  layer,  and  the  solution 
of  this  problem  is  determined  by  the  three  7r’s.  Definitions 
and  physical  interpretations  concerning  the  three  tv’s  are  as 
follows. 

7t\  indicates  a  generalized  resistance  to  oxygen  reduction, 
which  is  expressed  as  the  ratio  of  operating  cell  current 
density  to  exchange  current  density  (based  on  geometric 
area),  i.e. 

I  operating  cell  current  density 

7T[  =  -  ~  -  (17) 

( Avio)8c  exchange  current  density 


jtk  represents  a  normalized  resistance  to  proton  conduction, 
which  is  expressed  as  the  ratio  of  the  ohmic  voltage  drop 
across  the  catalyst  layer  to  the  Tafel-slope,  i.e. 

_  73c  _  I(8c/KeP)  _  ohmic  voltage  drop 

K  K^b  b  Tafel-slope 

and  7 To  is  a  normalized  resistance  to  oxygen  diffusion,  which 
takes  the  form  of  the  ratio  of  the  oxygen  consumption  rate 
to  the  diffusive  oxygen  transport  rate  at  a  concentration  gra¬ 
dient  of  Cq2/8c,  i.e. 


I8r 


7TD  = 


4 FC*02D% 


I/4F 


Do2(C*02/8c ) 


oxygen  consumption  rate 
diffusive  oxygen  transport  rate 


(19) 


In  order  to  understand  this  problem  better,  we  shall  not 
solve  Eqs.  ( 1 1)— ( 13)  directly.  The  coupled  system  of  three 
first-order  differential  equations  can  be  decoupled  by  differ¬ 
entiating  Eq.  (11)  with  respect  to  z,  and  substituting  di]/dz, 
dC/dz  into  the  so-obtained  equation  with  the  right-hand  side 
of  Eqs.  (12)  and  (13),  respectively.  This  then  becomes 

d 2i  f  i  —  1  4  d i 

5?  =  rD^+,,K'is 


(20) 
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Eq.  (20)  can  then  be  transformed  into  two  first-order  differ¬ 
ential  equations  by  the  introduction  of  an  additional  variable 
ip,  which  stands  for  the  dimensionless  reaction  rate.  We  now 
write 


df 

dz 


(21) 


dip 

dz 


jtd- 


C 


ip 


(22) 


Eqs.  (13),  (21)  and  (22)  represent  the  three  dependent  vari¬ 
ables:  C,  i,  and  ip.  Once  their  values  at  z  =  0  are  designated, 
the  nonlinear  first-order  initial-value  problem  can  be  solved 
to  yield  the  C,  i,  and  ip  profiles  using  numerical  methods 
(e.g.  the  Runge-Kutta  method).  We  already  have  the  initial 
values  of  i(0)  =  0  and  C(0)  =  1  from  Eqs.  (15)  and  (16). 
However,  the  value  of  ip(0)  must  be  chosen  to  ensure  that  the 


solution  profiles  satisfy  Eq.  (14).  Consequently  this  results 
in  a  shooting  problem  that  can  be  resolved  by  the  appropri¬ 
ate  shooting  technique  [16].  After  obtaining  the  solutions  to 
the  initial-value  problem,  the  profiles  of  C  and  ip  are  then 
substituted  into  Eq.  (11)  to  yield  the  spatial  variation  of  i). 


4.  Results  and  discussion 

It  is  necessary  to  remark  the  facts  that  the  dimension¬ 
less  parameter  n\  has  been  eliminated  through  the  deriva¬ 
tion  of  Eq.  (20),  which  was  then  transformed  into  Eqs.  (21) 
and  (22),  and  tti  is  not  involved  in  Eq.  (13).  Thus,  the  spa¬ 
tial  variation  of  the  dimensionless  reaction  rate  ip,  protonic 
current  density  i,  and  oxygen  concentration  C,  obtained  by 
solving  Eqs.  (13),  (21)  and  (22),  is  independent  of  tt\.  Their 
solution  profiles  are  determined  only  by  both  7Tk  and  tT[> 
This  interesting  phenomenon  signifies  that,  for  a  fixed  set 


Fig.  2.  Influence  of  7Tk  and  7Td  on  the  dimensionless  spatial  variation  of  (a)  reaction  rate,  (b)  current  density,  and  (c)  oxygen  concentration. 
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Table  1 


Electrode  properties  and  physical  parameters  for  base-case  conditions 


Electrode  property/physical  parameter 

Values 

Cell  temperature,  T  (K) 

353 

Cell  current  density,  /  (A/cm2) 

0.5 

Cathodic  transfer  coefficient,  ac 

2.0 

Reference  exchange  current  density  times  area. 

5  x  10~4  [1] 

Ayi0  (A/cm3) 

Catalyst  layer  thickness,  <5C  (cm) 

0.001 

Reference  oxygen  concentration,  CJ  (mol/cm3) 

4.62  x  10"6  [1] 

Effective  protonic  conductivity  in  catalyst 

0.002  [3] 

layer,  (S/cm) 

Effective  diffusion  coefficient  of  oxygen  in 

2x  10"4  [3] 

catalyst  layer,  (cm2/s) 

of  7Tk  and  ttu  values,  decreasing  the  value  of  tz\  (e.g.  the 
method  of  increasing  the  catalyst  density)  will  not  alter  the 
profiles  of  reaction  rate  /p,  protonic  current  density  i,  and 
oxygen  concentration  C,  but  fj  will  be  decreased  accord¬ 
ing  to  Eq.  (11).  Fig.  2  shows  the  influence  of  7Tk  and  7Td  on 
the  spatial  variation  of  ip,  i,  and  C  in  the  catalyst  layer.  The 
base  values  of  jti,  7Tk  and  7Td  are  1  x  106,  16.43,  and  1.402, 
respectively,  and  they  are  evaluated  from  the  base-case  con¬ 
ditions  given  in  Table  1 . 

It  can  be  seen  that  the  oxygen  concentration  profile  drops 
rapidly  (see  Fig.  2c)  as  the  value  of  7Td  is  increased  to  five 
times  its  base  value.  A  large  7Td  indicates  low  oxygen  per¬ 
meability,  so  the  oxygen  concentration  is  extremely  low  in 
the  right  portion  of  the  catalyst  layer.  Accordingly,  the  dis¬ 
tribution  of  the  reaction  rate  (see  Fig.  2a)  shifts  to  the  left, 
and  no  oxygen  reduction  occurs  in  the  right  portion.  As  the 
value  of  7Td  decreases  to  one-fifth  of  the  base  value  the  oxy¬ 
gen  will  penetrate  the  catalyst  layer  more  easily,  thus  the 
oxygen  concentration  profile  becomes  more  uniform  and  the 
distribution  of  the  reaction  rate  shifts  to  the  right. 

In  contrast  to  the  effect  of  changing  the  value  of  ttd,  in¬ 
creasing  the  value  of  jtk  to  five  times  its  base  value,  the  re¬ 
action  rate  in  the  left  portion  decreases  noticeably  because 
the  resistance  to  proton  conduction  is  increased.  Instead  of 
shifting  to  the  right-most  boundary,  the  reaction  rate  peaks 
around  z  —  0.8.  This  is  due  to  the  fact  that  the  oxygen 
diffusion  is  not  high  enough  to  provide  a  good  penetration 
through  the  catalyst  layer,  so  most  of  the  oxygen  is  con¬ 
sumed  somewhere  off  the  catalyst  layer/PEM  interface.  If 
we  reduce  the  value  of  jtk  to  one-fifth  of  the  base  value, 
the  resistance  to  proton  conduction  is  lowered,  hence,  the 
protons  are  transported  to  the  left  more  easily,  and  the  dis¬ 
tribution  of  reaction  rate  tends  to  level  off. 

Inserting  the  solution  profiles  of  C  and  ip  into  Eq.  (11) 
the  spatial  variation  of  the  dimensionless  overpotential  r/  is 
readily  obtained.  Fig.  3  gives  the  solution  profiles  for  the 
base  value  of  n\  associated  with  the  same  five  combinations 
of  7Tk  and  7Td  used  in  Fig.  2.  As  pointed  out  above,  when 
referring  to  Eq.  (6),  the  solution  profile  of  fj  is  a  negative, 
monotonically  decreasing  function.  Higher  values  of  7Tk  and 
7 td  indicate  greater  hindrance  to  the  transport  of  proton  and 


Fig.  3.  Spatial  variation  of  dimensionless  cathodic  overpotential  throughout 
the  catalyst  layer  for  the  same  combinations  of  tt^  and  ttd  used  in  Fig.  2. 

oxygen,  hence  more  voltage  losses  occur  across  the  catalyst 
layer.  As  can  be  seen,  for  a  given  value  of  tc\,  increasing  the 
values  of  7Tk  and  jt\-p  causes  a  greater  overpotential  across 
the  catalyst  layer,  and  vice  versa. 

The  main  purpose  of  this  study  is  to  investigate  the  influ¬ 
ence  of  7Ti,  7Tk  and  7Td  on  the  dimensionless  overpotential 
associated  with  the  whole  cathodic  catalyst  layer,  which  we 
designate  as  f;cat.  As  depicted  in  Fig.  3,  the  local  overpo¬ 
tential  varies  throughout  the  catalyst  layer.  In  this  study  we 
take  ?)ca t  =  r)(  1)  because  ij(  1)  is  equal  to  the  experimen¬ 
tally  measured  dimensionless  cathodic  overpotential  [3,6]. 
From  Eq.  (11)  and  the  foregoing  analysis,  it  is  perceived 
that  the  dimensionless  overpotential  associated  with  the  ca¬ 
thodic  catalyst  layer,  f)cat,  is  governed  by  jrj,  C(l),  and  fp(l). 
The  latter  two  are  in  turn  governed  by  7Tk  and  ttd.  We  de¬ 
duce  that  the  relationship  among  i)cat,  n\,  ttk  and  ttjj  can  be 
expressed  as 

Wat  =  -In  7t\  -  /(ttk,  ttd)  (23) 

The  influence  of  tc\  on  i)cat  can  be  decoupled  from  7Tk, 
jtd  and  takes  a  logarithmic  form,  whereas  the  influence  of 
both  7Tk  and  7Td  are  coupled  together.  The  first  term  on  the 
right-hand  side  of  Eq.  (23)  represents  the  activation  overpo¬ 
tential,  and  the  second  term  stands  for  the  mixture  of  ohmic 
overpotential  and  concentration  overpotential.  The  unknown 
function /(ttk,  n d)  can  be  determined  by  evaluating  i)cat  for 
a  given  value  of  tz\  associated  with  various  combinations  of 
ttk  and  ttd-  Fig.  4  shows  the  curves  of  f{n k,  7Td)  for  some 
constant  values  of  7To  with  respect  to  the  variation  of  7Tk. 
This  figure  presents  the  curves  of  tte>)  up  to  a  func¬ 

tion  value  of  60,  which  is  equivalent  to  0.913  V  in  the  mix¬ 
ture  of  ohmic  overpotential  and  concentration  overpotential 
(base-case  values,  ac  =  2.0,  and  T  =  353  K,  are  applied). 
This  range  is  wide  enough  to  cover  the  entire  operating  con¬ 
ditions  of  a  PEM  fuel  cell. 
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Fig.  4.  Mixture  of  ohmic  and  concentration  overpotential  as  a  function  of  tz k  and  .T|>. 

As  depicted  in  Fig.  4,  the  function  value  of  /(7Tk,  7Td) 
approaches  zero  as  both  jtk  and  : approach  zero.  Thus, 

Eq.  (23)  reduces  to  pcat  =  —In  n\  under  any  one  of 
the  following  conditions:  (i)  the  thickness  of  the  catalyst 
layer  is  extremely  thin,  (ii)  the  current  density  is  ex¬ 
tremely  low,  and  (iii)  both  the  effective  conductivity 
and  the  effective  diffusion  coefficient  are  extremely 
high. 

It  can  be  seen  from  Fig.  4  that  /(7Tk,  ttd)  is  an  increas¬ 
ing  function  of  jtk  and  jt [> .  In  Fig.  4,  one  can  see  that  the 
curves  for  7Td  >  20  coincide  with  that  of  jtd  =  20.  The 
reason  for  this  lies  in  the  limitation  of  resistance  to  oxygen 
diffusion.  With  the  value  of  7Td  approaching  20,  the  oxygen 
diffusion  becomes  very  slow,  and  the  oxygen  concentration 
profile  drops  abruptly  to  zero.  Thus,  the  reduction  of  oxy¬ 
gen  takes  place  only  within  an  extremely  thin  layer,  or  is 
even  restricted  to  the  left-most  boundary.  When  this  hap¬ 
pens,  increasing  the  value  of  i rjj  will  not  further  increase 
the  cathodic  overpotential.  On  the  other  hand,  the  curves  for 
ttd  <  0.02  coincide  with  that  of  7Td  =  0.02.  This  indicates 
that,  with  the  resistance  to  oxygen  diffusion  already  being 
very  small,  7Tk  is  the  only  controlling  factor  of  /(ttk,  7Td), 

Hence,  furthering  the  diffusion  capability  of  oxygen  will  not 
lower  the  cathodic  overpotential. 

An  empirical  correlation  can  be  developed  for  the  func¬ 
tional  form  of  /(7Tk,  ttd)  utilizing  the  curves  shown  in  Fig.  4. 

With  7Td  being  fixed  and  greater  than  1.2,  a  linear  relation¬ 
ship  between  the  function  value  and  jtk  is  observed.  Also, 
for  t td  being  smaller  than  0.8,  a  semi-logarithmic  relation¬ 
ship  is  found.  The  coefficients  of  correlation  can  be  obtained 


by  regression  analysis,  and  the  following  correlation  equa¬ 
tions  are  suggested: 

/  1  —  7TD  \ 

/<»K.»D)  =  (_01216_09(i35jJ>rK.  *d>1.2 

(24) 
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Fig.  6.  Fit  of  Eq.  (25)  to  the  calculated  function  values  for  ttd  <  0.8. 
and 

/(7Tk,  ttd)  =  (—0.65  —  0.37Td  +  2.2^d2)  +  1-025  In  7Tk, 

t rD  <  0.8  (25) 

Figs.  5  and  6  respectively  show  the  fit  of  the  correlation 
equations  to  the  calculated  function  values.  As  shown  in 
both  these  figures,  the  two  correlation  equations  are  in  good 
agreement  with  the  calculated  function  values.  The  maximal 
error  in  correlation,  which  occurs  in  the  low  7Tk  region,  is 
estimated  to  be  2.1  for  Eq.  (24),  and  0.65  for  Eq.  (25).  They 
are  respectively  equivalent  to  0.032  and  0.01  V  in  overpo¬ 
tential  when  the  base-case  values  of  cell  temperature  T  and 
transfer  coefficient  ac  are  applied. 

5.  Conclusions 

A  mathematical  model  for  the  catalyst  layer  of  a  PEM 
fuel  cell  cathode  has  been  developed,  the  formulated  equa¬ 
tions  were  nondimensionalized,  and  were  then  reduced  to 
a  system  of  dimensionless  differential  equations  associated 
with  three  dimensionless  parameters,  jt\,  jtk  and  7Td-  Here, 
n\  represents  the  resistance  to  oxygen  reduction,  7Tk  the  re¬ 
sistance  to  proton  conduction,  and  7Td  the  resistance  to  oxy¬ 


gen  diffusion,  hence  they  are  responsible  for  the  activation 
overpotential,  ohmic  overpotential,  and  concentration  over¬ 
potential,  respectively.  This  approach  has  simplified  the  pre¬ 
sentation  of  complicated  results  relevant  to  the  catalyst  layer, 
and  offered  a  useful  tool  for  the  investigation  into  the  PEM 
fuel  cell  cathode. 

The  profiles  of  oxygen  concentration,  reaction  rate,  and 
current  density  are  governed  only  by  7Tk  and  jtd-  They  influ¬ 
ence  the  cathodic  overpotential  through  the  spatial  variations 
of  reaction  rate,  and  oxygen  concentration.  The  influence  of 
7i]  on  )';cat  can  be  decoupled  from  7Tk,  jtd  and  takes  a  loga¬ 
rithmic  form,  whereas  the  influence  of  both  ttk  and  j rp  are 
coupled  together.  The  dependence  of  the  mixture  of  ohmic 
overpotential  and  concentration  overpotential  on  ttk  and  jtd 
was  depicted  in  Fig.  4,  which  can  be  conveniently  referred 
to  when  assessing  the  performance  of  a  PEM  fuel  cell. 

Further  more,  an  empirical  correlation  was  developed  us¬ 
ing  the  results  presented  in  Fig.  4.  This  has  yielded  the  fol¬ 
lowing  two  correlation  equations: 

„  /  1  —  7TJ)  \ 

'Jc.,  —  —In  +  (o.1216  +  0.9635jtd /  JIKl  ”D-h2 

(26) 

i)cat  =  —In  jti  +  (0.65  +  0.37Td  —  2.27Td2)  —  1.025  In  ttk, 

ttd  <  0.8  (27) 

This  correlation  is  considered  to  be  helpful  for  the  analysis, 
prediction  and  optimization  of  the  performance  of  a  PEM 
fuel  cell. 
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